Deterministic Generation of Arbitrary Symmetric States and Entanglement Classes 
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We propose a method to generate arbitrary symmetric states of N qubits, which can be easily 
associated with their entanglement classes. It is particularly suited to quantum optics systems like 
trapped ions or superconducting circuits. We encode each qubit in two metastable levels of the 
system and use a bosonic quantum bus for creating the states. The method is deterministic and 
relies on a sequence of selective unitary gates upon the qubits within the system coherence time. 
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Entanglement is one of the most significant features 
in quantum information science. It is considered crucial 
for achieving an exponential speedup in quantum com- 
puting protocols and an essential resource in quantum 
communication pQ. Even though a full characterization 
of entanglement is still lacking [2j [3] , several important 
steps have been achieved in order to classify entangled 
states according to different criteria. Among these, the 
most prominent is the classification under Stochastic Lo- 
cal Operations and Classical Communication (SLOCC). 
This was firstly obtained for the 3-qubit case [4 and 
subsequently extended to the 4-qubit case [5-8 and 
larger numbers of qubits [9 . Two 7V-qubit pure states 
|^) and |0) are equivalent under SLOCC if and only if 
|^) = Ai®A2®...®Ajsf\(j))^ where {Ai} i= i^^N are invert- 
ible matrices, not necessarily unitary. Two states belong 
to the same SLOCC entanglement class if and only if they 
are equivalent under SLOCC. The motivation for this en- 
tanglement classification is practical, i.e., two states be- 
longing to the same SLOCC class are able to perform the 
same quantum information tasks up to a finite probabil- 
ity. Extensions to N > 4 qubits are cumbersome, due to 
the fact that the number of entanglement classes becomes 
infinite already for N = 4 [4 . On the other hand, it is 
sensible to group the infinite number of classes in a cer- 
tain set of families [SJ [10] , containing classes that share 
similar mathematical properties. It is also interesting to 
analyze entanglement classes and families of states that 
have specific symmetries. In this respect, the classifica- 
tion for 7V-qubit pure states that are symmetric under 
qubit permutation has been achieved jT0HT2] . We point 
out that this set of symmetric states contains an infinite 
number of SLOCC classes in some of the entanglement 
families. The experimental access to the vast majority of 
states in the symmetric classes remains an open problem, 
whose solution would enhance the use of multipartite en- 
tanglement in quantum information protocols. 



On the other hand, an increasing interest for opera- 
tional entanglement classifications, with direct connec- 
tion to experimental setups, has given rise to propos- 
als based on atoms or ions coupled by spontaneously 
emitted photons and interference in the detection pro- 
cess [13 . However, this kind of proposals have a small 
success probability that decreases exponentially with the 
number of qubits, restricting their scalable implementa- 
tion [14 . Trapped ions are one of the most controllable 
quantum systems [T5l [16] and, nevertheless, a limited 
set of prototypical entangled states such as W [17] or 
Greenberger-Horne-Zeilinger (GHZ) [18] states of up to 
14 qubits have been achieved [19]. Proposals for symmet- 
ric Dicke state generation, but not their superpositions as 
required for the much larger set of arbitrary symmetric 
states, have also been put forward [20-26 . An alter- 
native quantum platform, as superconducting quantum 
circuits, is the most promising solid-state quantum tech- 
nology [27[ |28] . Significant steps towards the generation 
of multipartite entangled states in this implementation 
have been given [29H32] . However, a protocol for gener- 
ating the broad set of arbitrary 7V-qubit symmetric states 
in these or other technologies is still missing. 

In this Letter, we propose an efficient protocol for gen- 
erating arbitrary symmetric 7V-qubit pure states in quan- 
tum optics systems like trapped ions or superconducting 
circuits. Our protocol is deterministic and feasible with 
current technology, making use of selective interactions 
to map product states with fixed number of excitations 
to symmetric Dicke states with the same number of ex- 
citations. We also propose a mapping between symmet- 
ric entanglement classes and experimental parameters in 
quantum controllable systems. 

We first introduce the basic idea of our protocol and 
subsequently discuss implementation details. We con- 
sider a set of N qubits that will be coupled with a bosonic 
mode acting as a quantum bus. The latter can be a mo- 
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tional mode in the trapped-ion case, or a microwave field 
in the superconducting circuit case. Our aim is to gener- 
ate an arbitrary symmetric state of the N qubits, which 
can be expressed as an arbitrary superposition of sym- 
metric Dicke states with variable number of excitations, 



\^n)s 
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k=0 



Ck\DN,k): 



(i) 



where \Dn^) is the symmetric Dicke state with N qubits 
and k excitations (we will give a 4-qubit example below). 
We point out that one can directly map the Ck coefficients 
in this state to the corresponding symmetric entangle- 
ment classes through Vieta's formulas [TQl [TTJ [13] . We 
will use this property at the end of the paper to establish 
a correspondence between entanglement classes for sym- 
metric states and experimental parameters in quantum 
optics systems. 

As the starting point of the protocol, we prepare the 
N qubits in states \g). Then, we generate the same su- 
perposition as in Eq. ([!]) , but in the basis of states having 
k last qubits in state \e) and the remaining first N — k 
qubits in state \g), 



N 



(2) 



k=0 



where \(k)) = \gig2---9N-k^N-k+i--^N) • To illustrate 
how this task could be implemented let us consider a 
particular case of a 4-qubit state. We start from the 
state with zero excitations |^o) = \d99d)- Let us apply 
a local rotation on the fourth qubit such that we have 
= a \gggg) + Po\ggge). Then, we apply a rotation 
on the third qubit, conditioned to the fourth qubit being 
in the excited state |e). The way to do this in terms of 
elementary gates is shown in Ref. pQ. This leads to state 
1^2) = ao\gggg) + Pooti\ggge) + p Pi\ggee). Implement- 
ing a similar operation on the second qubit conditioned 
to the state of the third qubit, and a rotation on the first 
qubit conditioned to the state of the second, we finally 
obtain the desired four qubit state, 

1^4) = a \gggg) + Poa^ggge) + p Pia 2 \ggee) 

+ PoPiP2Cts\geee)+p p 1 p 2 Ps\eeee). (3) 

Coefficients and fij from Eq. ([3| can always be found 
out of coefficients c/~ in Eq. Following this procedure, 
we can generate an arbitrary superposition of states with 
k excited qubits as defined above. The needed elemen- 
tary gates grow polynomially in the number of qubits of 
the target symmetric state. With this example, we have 
given a specific protocol for generating ^4). 

The next step in the protocol is to coherently map 
the states with k excitations to symmetric Dicke states 
\DN,k) of N qubits with k excitations. As a first strategy 
to undertake this task, we observe that a symmetric Dicke 
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FIG. 1: (color online), (a) Scheme of multiqubit system as 
proposed in the text. In one step of the proposed protocol, 
the first n qubits are addressed with a global red sideband 
interaction of Rabi frequency g\ and detuning Ai, while the 
last qubit (first one from the right in the figure) is addressed 
with a second red sideband driving of Rabi frequency #2 and 
detuning A2. Here |e, 0) (resp., \g, 1)) is the excited internal 
state with no excitations in the bosonic mode (resp., ground 
internal state with one excitation in the bosonic mode), (b- 
d) Probability amplitudes for the evolved state of being in 
the corresponding states shown in the curves, as a function 
of time and the respective step, in the protocol for obtaining 
state |^ 4 )s in Eq. (l), with c k = {1, 2, 3, 4, 5}/^55. Here 
ti=tt/(4/3). 



state with k excitations can always be decomposed as 
follows, 



dN,k 

d>N-l,k 

d>N,k 



D N - lik )\g), 



(4) 



where dw,k is the normalization of a symmetric Dicke 
state of N particles with k excitations, e.g., \Dn,i) = 
(\eg...g) + \geg...g) + ... + \g...ge))/d N ^. From this ex- 
pression, we learn that a symmetric Dicke state of N 
particles with k excitations can be obtained by symmet- 
rically adding one excitation to a symmetric Dicke state 
of N — 1 particles with k — 1 excitations. Thus, if we are 
allowed to realize a Rabi flopping between |£)jv-i,fc-i)| e ) 
and \DN-i t k)\g), we can obtain a Dicke state of N par- 
ticles with k excitations. Such Rabi flopping could be 
realized by an interaction that annihilates one excitation 
in one qubit and distributes it symmetrically in a set 
of qubits. This can be achieved by coupling a set of n 
qubits with a bosonic quantum bus by means of a red- 
sideband interaction with a detuning Ai and a collective 



3 



Rabi frequency g\. Additionally, we need another qubit 
that is coupled with a second red-sideband driving with 
a detuning A2 and Rabi frequency #2, see Fig. Pro- 
vided that Ai,A2 ^> QxiQi-, the system is described by 
the Hamiltonian 



i^j i=l 



St 



P*ajae- M ), (5) 



where G{ is the spin ladder operator associated to qubit 
z, and a is the respective ladder operator associated to 
the additional qubit, 5 = Ai — A2 — A2 + Ai, f3 = 
g$g 2 /A and A, = |ft| 2 /A„ with 1/A = (1/2)(1/Ai + 
I/A2). The main advantage of this Hamiltonian is 
that it allows independent evolution inside the sub- 
spaces {\D n ^-i)\e), \D nj k)\g)}, and it permits to selec- 
tively choose which subspace will evolve, provided by an 
appropriate tuning of the parameter 5. The selectivity of 
a given subspace is possible when Ai ^> j3. The mapping 
of the 7V-qubit basis, as defined in state ([2|, can be un- 
derstood from our example for the four qubit state (|3|. 
As shown in the following diagram, 
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we access the first two qubits in such a way that, un- 
der the dynamics given by Hamiltonian ([5|, the system 
only evolves within the subspace {\Di$eee), \Di^gee)}. 
Thus, after an evolution time such that the corresponding 
amplitudes in Q are satisfied, we obtain state \D 2 ^ee). 
In steps (ii) and (iii), S is chosen in such a way that 
only \D2,iee) or \D2^ee) will respectively evolve to states 
1^3,2^) and \Ds^e). In Table [l| the different parameters 
for each step are shown. We also include the fidelities 
for the respective steps obtained performing numerical 
simulations of the dynamics, for obtaining state ^4)5 in 
Eq. ([l]), with Ck = {1, 2, 3, 4, 5}/VbE. Assuming that the 
fidelity of obtaining the initial state ([3| is close to 1, the 
total fidelity for generating state ^4)5 is ~ 0.989. The 
fidelity loss is due to the dispersive character of the in- 
teraction in Eq. ^ and imperfections in the selectivity 
condition. The probability amplitudes of the dynamics 
of the process can be seen in Figs, [I}}, lc, and Id. 

Let us now give details of possible implementations. 
Any physical system that allows for detuned red-sideband 
interactions among well characterized qubits is appro- 
priate for our proposal. Some examples are trapped 
ions [I5j [16] , or superconducting qubits [27j [28] . We con- 
sider now a possible trapped ion implementation. The 
scheme could be implemented, e.g., using a linear string 
of ions in a Paul trap. To generate a symmetric 7V-qubit 
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TABLE I: Parameters for the preparation of state |^4)s in 
Eq. {l}, with c k = {1, 2, 3, 4, 5}/VtE, in units of g 1 . 



entangled state N ions are loaded into the trap. Stan- 
dard techniques of Doppler and resolved sideband cool- 
ing prepare a linear ion crystal, with the axial center-of- 
mass (COM) motional mode cooled close to the ground 
state. The method described to prepare the initial state 
in Eq. and its generalisation, can be done using the 
universal trapped- ion gate set in Refs. [33 , 34 . To derive 
an estimate of the four-ion protocol time, we consider 
singularly ionised calcium ( 40 Ca + ) ions, which are used 
by many groups around the world. The ions provide a 
long-lived (~ 1 s) quadrupole transition at optical fre- 
quencies (729 nm), between the |5i/ 2 ) and \D 5 / 2 ) Zee- 
man states. A weak static magnetic field is applied to 
lift the Zeeman degeneracy, providing a large number of 
transitions that can be used to encode quantum infor- 
mation. With typical magnetic fields of a few Gauss, the 
splitting is around 7 MHz between the ^5/2) levels and 
10 MHz between the IS1/2) levels. For a four ion string, 
Rabi frequencies can typically be tuned anywhere from 
to about 1 MHz, and red-sideband interactions using 
the center-of-mass axial mode are straightforwardly im- 
plementable. Using g\=2ir x 20 kHz we calculate a time 
of 5 ms to prepare the N = 4 ion Dicke state in Eq. 
with Ck = {1, 2, 3, 4, 5}/ a/55- Current trapped- ion tech- 
nology achieves hundreds of gates per experiment inside 
a total coherence time of about 30 ms [34]. Our proto- 
col will allow to generate an arbitrary N = 4 symmetric 
state in trapped ions with fidelities of about 0.99 inside 
this coherence time. 

In order to detect the final state ^4)5 after the pro- 
tocol, one may apply recently developed techniques for 
symmetric state tomography [35] . The advantage of these 
techniques is that the number of observables one has to 
measure to perform the full tomography grows only poly- 
nomially in the number N of qubits. This is due to the 
symmetry of the considered states under permutation of 
the qubits. To perform this tomography it is enough to 
realize a polynomial number of local rotations and cr z 
measurements upon the different qubits [35]. The latter 
can be done in the trapped-ion example via resonance 
fluorescence with a cycling transition to a metastable ex- 
cited state [15] [16] . 

We remark that this is a deterministic protocol, i.e., it 
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happens with probability 1, up to experimental imper- 
fections. This should be compared with existing proba- 
bilistic proposals and experiments for generating entan- 
glement via photon emission and detection, with success 
probability for a two- ion entangled state of ~ 10 -8 [T4] . 
The success probability for N = 4 ions is, indeed, much 
smaller in these probabilistic methods given that it de- 
creases exponentially with N. 

Another appeal of our proposal is the mapping be- 
tween symmetric entanglement classes and experimental 
parameters in quantum optics systems. Given that this 
protocol enables the generation of arbitrary symmetric 
states of TV qubits, we can in particular generate all sym- 
metric entanglement classes [10]. One may establish a 
link between each entanglement class and the experimen- 
tal configuration of the specific implementation. First 
of all, one establishes a correspondence between the Ck 
coefficients of a specific I^Ar)^ in Eq. ([!]) and the ex- 
perimental configuration to generate it. Indeed, the 
coefficients are already determined by the set of gates 
necessary to generate state \^n) hi Eq. In our spe- 
cific example, the set of parameters c^, i = 0, 1, 2, 3, 
associated with local gates and conditional rotations pQ, 
determines the numbers, and it is attainable, in the 
trapped-ion example, with the gate toolbox introduced 
in [33] [34]. Then, one maps the Ck coefficients to the 
corresponding entanglement class through Vieta's formu- 
las nni nu us] . 

Instead of encoding the initial superposition Q of 
states with different number of excitations with the set of 
ket vectors one could encode it on the set of Fock 

states associated with the bosonic mode, \k) pn . Thus, 
this other complementary approach requires the initial- 
ization of an arbitrary bosonic state with up to N excita- 
tions. In Ref. [25 , it was shown how to create arbitrary 
Dicke states in trapped ions. The idea there was to first 
generate a /c-phonon motional Fock state of, e.g., the 
center of mass mode, |fc) pn , and subsequently to transfer 
it via adiabatic passage with a red sideband transition 
to the symmetric Dicke state with k excitations, \Djy^)- 
Notice that in this way only Dicke states, and not ar- 
bitrary symmetric states, are generated. To achieve the 
latter, we suggest to first initialize the center of mass mo- 
tional state in an arbitrary superposition of Fock states 
with up to TV excitations, ^2 k Ck\k) pn . In order to attain 
this, we refer to the proposal by Law and Eberly [36] 
for creating an arbitrary electromagnetic field in cavity 
QED using Jaynes-Cummings [37] and carrier interac- 
tions. This proposal is directly applicable to trapped 
ions by substituting the electromagnetic field with the 
motional center of mass mode. Finally, we propose to 
apply an adiabatic passage with the red sideband tran- 
sition, which will map each phonon Fock state with k 
excitations to each symmetric Dicke state with k excita- 
tions. This will produce the desired arbitrary superposi- 
tion of symmetric Dicke states, ^2 k Ck\DN,k), i- e -? an ar- 



bitrary symmetric state of N qubits. The same idea can 
be applied to superconducting qubits by substituting the 
arbitrary motional state by an arbitrary microwave-field 
Fock state, which can be, in principle, similarly achieved. 

Summarizing, we have introduced a deterministic pro- 
tocol for generating arbitrary symmetric states of TV 
qubits, and their respective entanglement classes, by 
means of selective interactions. This will allow the cre- 
ation of this kind of states for a few ions or superconduct- 
ing qubits with high fidelities. We have also given a com- 
plementary approach based on the initialization of an ar- 
bitrary superposition of Fock states and subsequent map- 
ping to the equivalent Dicke state superposition. These 
proposals will contribute to enhance the sets of states 
that can be achieved in controllable quantum systems 
with current technology. 
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